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Abstract 

We present a pedagogical overview of recent theoretical work on unconventional quantum phases and quantum 
phase transitions in condensed matter systems. Strong correlations between electrons can lead to a breakdown 
of two traditional paradigms of solid state physics: Landau's theories of Fermi liquids and phase transitions. We 
discuss two resulting "exotic" states of matter: topological and critical spin liquids. These two quantum phases do 
not display any long-range order even at zero temperature. In each case, we show how a gauge theory description 
is useful to describe the new concepts of topological order, fractionalization and deconfinement of excitations 
which can be present in such spin liquids. We make brief connections, when possible, to experiments in which the 
corresponding physics can be probed. Finally, we review recent work on deconfined quantum critical points. The 
tone of these lecture notes is expository: focus is on gaining a physical picture and understanding, with technical 
details kept to a minimum. 
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1. Introduction 

Traditional condensed matter physics is shaped 
by two ideas that can be traced back to Landau: the 
Fermi liquid theory and the theory of phase tran- 
sitions [1] . These principles have proven very fruit- 
ful and enabled great progress in describing simple 
metals, semiconductors and insulators. However 
starting from the experimental discovery of the 
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fractional quantum Hall effect [2] , and before that 
theoretical predictions in Id quantum spin chains 
[3] , it has become clear that a large number of phe- 
nomena observed in condensed matter materials 
lie beyond the description of these two fundamen- 
tal paradigms. Such effects are observed in systems 
with strongly correlated electrons, which include 
for instance high T c superconductors [4], heavy 
fermion materials [5] or transition metal oxides [6] . 
In these systems the interactions between the elec- 
trons need to be considered explicitly and are pri- 
marily responsible for the observed phases of mat- 
ter, which cannot be explained in a Fermi liquid 
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picture. These new quantum mechanical phases are 
present at zero temperature, but their effects can 
also be felt at higher temperatures. Furthermore 
the phases of matter which emerge do not neces- 
sarily differ by symmetry, although quantum phase 
transitions between them are sometimes possible: 
this is not reconciliablc within the Landau theory 
of phase transitions. Therefore the standard Lan- 
dau paradigms do not apply and one must seek 
new approaches to understand the nature of these 
exotic phases. 

These notes are based on the lectures entitled 
"Exotic quantum phases and phase transitions in 
correlated matter" given by one of us (M.P.A.F) at 
the FPSPXI summer school in Leuven. They are 
meant to be pedagogical and they also reflect the 
process of understanding of the scribes (F.A. and 
A.M.W.). The lecture notes can therefore be useful 
as a basic introduction to the topic. We have also 
tried to keep in the tone of these notes the lively 
spirit of the lectures. 

The discussion will consider one of the sim- 
plest systems in which strongly correlated electron 
physics emerges, that of Mott insulators and will 
focus in particular on the example of spin liquid 
phases. Spin liquids are exotic phases that occur in 
quantum spin systems with odd number of spins 
S = 1/2 per unit cell when no symmetries are 
broken at zero temperature; as such, they can be 
thought of as unusual quantum paramagnets. We 
start with a discussion in Sec. 2 of the situations 
where the interactions between electrons matter 
and how these strong correlations can break down 
the traditional solid state physics paradigms. In 
the following two sections, we present two known 
classes of systems where these strong interactions 
stabilize exotic phases: topological (Sec. 3) and 
critical or algebraic (Sec. 4) spin liquids. We then 
present in Sec. 5 recent work on an explicit exam- 
ple of a quantum phase transition which cannot be 
explained naturally in terms of a Landau theory, 
and wc conclude in Sec. 6. 



2. Breakdown of traditional paradigms 

2.1. Fermi liquid theory and strong correlations 

Fermi liquid theory has been remarkably suc- 
cessful explaining why band theory works, even in 
systems with interacting electrons. Band theory 
basically assumes free electrons. Assuming a one- 
to-one correspondence between the states of the 
free system and those of the interacting one, the 
Fermi liquid theory shows that interactions essen- 
tially dress up the electrons on the Fermi surface 
into a "quasi-particle" with admittedly renormal- 
ized characteristics (effective mass etc.), but with 
retained independent properties - justifying band 
theory. 

According to the resulting band theory, elec- 
trons have allowed energy levels (bands) separated 
by forbidden energy regions (gaps) in the presence 
of the ionic periodic potential. For temperatures 
smaller than needed to thermally excite an electron 
over the gap, the electrons cannot move around if 
the highest energy band is filled and consequently 
the material is an insulator, whereas if the highest 
band is partially filled the material is a metal. Ef- 
fectively, band theory predicts that systems with 
an even number of electrons per unit cell are usu- 
ally insulators and those with an odd number are 
always metals. 

These predictions are correct if the bands are 
wide enough, as compared to the characteristic 
Coulomb interaction between electrons, such that 
the Fermi liquid assumption holds. This is usu- 
ally the case for Fermi surfaces formed by electrons 
with an atomic s or p character. For materials with 
electrons at the Fermi energy from partially filled 
atomic d or / orbitals, the overlap of these orbitals 
on neighbouring atoms is usually much smaller, 
since the interatomic spacing is still determined by 
the s and p orbitals. This results in a smaller kinetic 
energy of hopping and therefore narrower bands. 
In such cases the electrons are said to be strongly 
correlated. It is worth stressing that the Coulomb 
energy has the same magnitude in either case - it is 
the kinetic energy of the electrons which differs be- 
tween systems of strongly correlated electrons and 
those correctly described by Fermi liquid theory. 
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The validity of Fermi liquid theory is determined 
by the ratio of Coulomb and kinetic energies. 

Strong interactions lead to a correlated motion 
of electrons, and understanding its precise nature 
is the challenge faced by theories of strongly cor- 
related electron systems. For larger magnitudes of 
Coulomb repulsion, the electrons arc less likely to 
hop to a neighbouring site and they effectively self- 
localize on a given lattice site. To be more explicit, 
the electrons can virtually hop to another site and 
back, however the materials are insulators and will 
not conduct. Such materials, with an odd number 
of electrons per unit cell, that do no conduct be- 
cause of Coulomb repulsion arc named Mott insu- 
lators, as opposed to band insulators which result 
from filled energy bands. In fact band theory, which 
ignores the correlation effects that determine the 
nature of these materials, predicts that they are 
metals. Examples of such materials are transition 
metal oxides (e.g. cuprates, manganites, chlorides) 
with partially filled 3d and Ad bands and rare earth 
and heavy fermion materials with partially filled 
4/ and 5/ bands. 

The simplest model which takes into account the 
correlation effects is the Hubbard model 

H = -tJ2 [ c \, a c i,<* + h.c.) + U^n iA n i , i , (1) 

(ij).a i 

in which the inter-site hopping of the electrons 
t competes with on-site repulsion U. Here cj, Cj, rij 
are respectively creation, annihilation and density 
operators of an electron at lattice site i with spin 
index a which can take on up or down values (a =| 
, I). The model considers one electron per cell on a 
lattice and (ij) denotes pairs of neighbouring sites 
i and j. Ift/U is large, electrons are essentially free. 
In the other limit, U/t ^> 1, the electrons are site 
localized and the nature of the system is shaped 
by the residual spin physics. 

Kinetic energy favours two neighboring site- 
localized electrons to have antialigned spins, as 
a hop to a neighbouring site with the same spin 
orientation is not allowed due to the Pauli exclu- 
sion principle. Therefore at half-filling and in the 
large repulsion limit, the most elementary model 
that captures the dominant interactions in Mott 
insulators is the S = 1/2 Heisenberg Hamiltonian 
with antiferromagnetic exchange J 



if = J]T S ; - Sj. (2) 

(«> 

Here S = (S x ,S y ,S z ) are S = 1/2 operators, 
which (up to a constant) are equal to the Pauli 
matrices (a x ,a y , <r z ). The strength of the exchange 
constant is determined in terms of the original en- 
ergy scales by a second order virtual hopping pro- 
cess: J ~ t 2 /U. 

Although these models are quite simple to write 
down, they prove very complicated to solve. The 
Hubbard model, as well as a majority of the quan- 
tum spin models, may not be solved analytically 
apart from a few special cases. Moreover, unlike 
in classical problems, numerical calculations are 
generally not helpful. The exponentially large size 
of the Hilbcrt space prevents an exact diagonal- 
ization of H and the fermion sign problem renders 
Monte Carlo calculations inefficient. Since the 
fermion wavefunction changes sign under inter- 
change of particles, calculating the partition func- 
tion requires summing a large number of positive 
and negative terms, which results in large numeri- 
cal error. Therefore to solve the problem one must 
try to deduce behaviours that might appear based 
on toy models and/or use basic physical notions 
such as universality. In classical systems, univer- 
sal, system independent physics emerges when the 
correlation length becomes much larger than the 
interatomic spacing. The difference in the order of 
the magnitudes of the two characteristic lengths 
is responsible for a length scale separation and 
leads to system independent behaviour. In quan- 
tum system, the scale separation necessary for the 
existence of universal phenomena, is provided by 
the large characteristic energy scale of electrons in 
solids (t or U in the Hubbard model), which is or- 
ders of magnitude larger than room temperature. 

2.2. Ordering in a Mott insulator 

Let us now focus on the behaviour at zero 
temperature of Mott insulators, which can be de- 
scribed by a spin Hamiltonian such as the Heisen- 
berg model Eq. (2) and its variants. 

At T = 0, Mott insulating systems usually break 
symmetries and develop order. For exemple, the 2d 
Heisenberg model on the square lattice develops 
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Fig. 1. The antiferromagnetic (or Neel) state on a square 
lattice. 
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Fig. 2. (a) A valence bond (or a dimer) is a singlet pair 
formed by two neighbouring spins on a lattice, (b) The 
Valence Bond Solid (VBS) formed by crystalline ordering 
of these dimers. 

antiferromagnetic (Neel) long range order associ- 
ated with the breaking of spin rotation symmetry 
(see Fig. 1). Another possibility is that two neigh- 
bouring spins prefer to first pair themselves in a 
singlet ^7|(|TI) — UT})> forming a valence bond (a 
"dimer"). This is denoted in Fig. 2A by an oval. 
These dimers can then order amongst themselves, 
forming a valence bond solid (see Fig. 2B) which 
breaks translational symmetry. In both cases, these 
phases can be described by a local order parame- 
ter. Transitions in or out of these phases are then 
formulated within Landau's traditional theory of 
phase transitions, which expands the free energy 
in powers of this order parameter. 

Our interest in the following sections lies else- 
where. We want to consider 5 = 1/2 systems which 
do not break any symmetries when cooled to zero 
temperature - they do not possess long-range or- 
der. These systems, dubbed spin liquids, possess 
however some kind of exotic order, which cannot 
be understood in terms of a local order param- 



eter. A theorem recently proven by Hastings [7] 
guarantees the existence of exotic quantum ground 
states in Mott insulators with no broken symme- 
tries. This theorem, a generalization of the Lieb- 
Schultz-Mattis theorem [8] to dimensions larger 
than 1, states that for 5 = 1/2 systems with one 
spin per unit cell on a two-dimensional lattice with 
periodic boundary conditions (on a torus), in the 
absence of symmetry breaking the ground state is 
separated from the first excited state by an en- 
ergy gap that vanishes in the thermodynamic limit: 
Ei - E < for a system of size L by L. This 

means that in the thermodynamic limit, there can- 
not be a singly-degenerate quantum paramagnet 
with a finite energy gap (see Fig. 3a). In the follow- 
ing sections, we will describe two ways out of this 
situation and the exotic properties of the associ- 
ated systems. Firstly we discuss in Sec. 3 topologi- 
cal spin liquids, where the disordered ground state 
is degenerated (see Fig. 3b) according to the topol- 
ogy of the system. These systems have an emer- 
gent gauge structure and admit excitations with 
fractional quantum numbers. In Sec. 4, we will de- 
scribe critical or algebraic spin liquids, with gap- 
less excitations above a ground-state that possesses 
power-law correlations (see Fig. 3c) . Such systems 
behave without fine-tuning as if they were at a crit- 
ical point. 



3. Topological spin liquids 

3.1. Description of topological spin liquids 

As their name suggests, topologically ordered 
spin liquids have peculiar behaviour associated 
with the topology of the system. They have a de- 
generate ground state in the thermodynamic limit 
with a degeneracy that depends on the topology 
of the system. The ground state wavef unctions, 
associated with the degenerate states, cannot 
be distinguished by any local measurement but 
are globally distinct. These states also have an 
emergent gauge structure which supports gapped 
particle-like excitations, with quantum numbers 
that cannot be built up from electron quantum 
numbers. It is possible to partially classify such 
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Fig. 3. Possible energy spectra (and associated real-space cartoon picture) of a spin liquid state with a disordered ground 
state in the thermodynamic limit, (a) For 5=1/2 spin systems with odd number of spins per unit cell, the Hastings' theorem 
forbids non-degenerate ground state with a finite energy gap. (b) Topological spin liquids possess a multiply-degenerate 
ground-state (with a topology-dependent degeneracy) with a gap to all excitations in the bulk. The cartoon picture presents 
a RVB liquid state of dimers with a delocalized free spinon S = 1/2 excitation. See Sec. 3. (c) Critical spin liquids possess a 
ground-state with power-law correlations and gapless excitations. In real space, critical spin liquids can typically be seen as 
valence bonds on many length scales. The excitations are spinons interacting with this valence-bond background. See Sec. 4. 



spin liquid phases in terms of the symmetry group 
of the emergent gauge theory. General discussions 
on topological order and its classification can be 
found in the book of Wen [21]. For simplicity, we 
will focus here on the simplest example of a topo- 
logically ordered system: a spin liquid with Zi 
topological order on a 2D square lattice. 

As proposed by P.W. Anderson [9], spin liquids 
can appear as a result of melting a valence bond 
solid. Consider a generalized Heisenberg model 
with not only nearest neighbours interactions 
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with Jij > 0. As already mentioned, in some cases, 
two neighbouring spins will form, due to the anti- 
ferromagnetic interaction, a singlet pair (a valence 
bond). These singlets can possibly align to form 
the valence bond solid (VBS), depicted in Fig. 2b. 
The bonds between four neighbouring sites can also 
fluctuate quantum mechanically between two pos- 
sible positions, as shown in Fig. 4a. These reso- 
nances between the bonds can lead to the melting 



Fig. 4. (a) The two resonating states connected by quantum 
fluctuations on a plaquette. The flux operator F p flips the 
bonds between the two resonating states (see Sec. 3.2). (b) 
A Resonating Valence Bond (RVB) liquid is a disordered 
state of valence bonds. 

of the solid to another state of matter with no long- 
range correlations of any kind: a resonating valence 
bond (RVB) liquid (Fig. 4b). The RVB liquid can 
be described using a wavefunction which is a su- 
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perposition of many valence bond configurations. 

What are the excitations of these two different 
phases? A valence bond may be broken, at an en- 
ergy cost of the order of J, to create an S = 1 
(gapped) excitation in which the spins are aligned 
(see Fig. 5a). Let's try to separate these two S = 
1/2 "spinons". The energy cost for this separa- 
tion will grow linearly with separation in the VBS, 
because there will be a mismatch in dimer alig- 
ment all along the path separating the two spinons 
(see Fig. 5b). The spinons are therefore confined 
in the VBS, similarly to quarks in elementary par- 
ticles. They cannot exist as finite energy excita- 
tions. However in the RVB liquid state, the re- 
gion between the two spinons is not a defect line, 
as in the solid, therefore the energy costs stays fi- 
nite with separation. The reorganization of the va- 
lence bonds to accomodate the extra spinon is lo- 
cal - or at least the energy cost is local. Therefore 
the spinons are deconfined in the liquid state (see 
the cartoon picture of Fig. 3b). Furthermore the 
spinons carry the S — 1/2 spin of the electron, but 
do not carry any additional electrical charge. The 
charge density in the topological spin liquid is uni- 
form and equal to one electron per cell. The exis- 
tence of these fundamental excitations which can- 
not be described by electron quantum numbers is 
referred to spin-charge separation. Quantum num- 
bers are said to be fractionalized, since the spin 
S = 1/2 of the spinon is half the one of the con- 
ventional spin- flip excitation with S = 1. 

3.2. Gauge theory formulation 

Due to their physical origin, the spin interactions 
in the Hamiltonian (3) are typically local. There is 
no explicit long-range interaction. How can we then 
understand the confincmcnt/dcconfmcment prop- 
erties of the spinons, which might be far apart? 
The natural framework for this is gauge theory. Let 
us now formulate the quantum dynamics of the va- 
lence bonds in such terms. 

To formulate the problem in terms of a gauge 
theory, consider a "spin" , on each lattice link, be- 
tween two sites i and j. This "spin" is described 
by a set of Pauli matrices er^ , the x component of 
which takes on the value erf.- = 1 if there is no va- 




Fig. 5. (a) Breaking of a valence bond results in two spinon 
excitations, (b) In a VBS, two separated spinons have cre- 
ated a mismatch in the arrangement of dimers: the energy 
cost will be linear in the distance between spinons. 

lence bond (no dimer) on the given link and 
crfj = — 1 if there is one. We want to define an op- 
erator which, when acting on four spins on a pla- 
quette, flips the bonds between the two resonating 
states, as depicted in Fig. 4a. This is easily done 
using operators that anticommute with each other, 
such as Pauli matrices. The bond flipping opera- 
tor, called the plaquette flux operator, is defined on 
the four sites k, I around the plaquette as F p — 
o~ij°~jk a ki ~ii- The Hamiltonian which describes the 
resonating valence bonds is a sum over all plaque- 
ttes of the flux operators H = —KJ2 p (ijki) Up 
to this point we have not taken into account the 
constraint that each spin can form only one sin- 
glet, i.e. only one valence bond comes from each 
site. Such a local constraint on the Hilbert space of 
this "quantum dimer model" may be encoded in a 
gauge charge Qi for each site i. The gauge charges 
are defined as the product of the bonds at a given 
site and we constrain them to be equal to — 1 , Qi = 
Ilj=i a ij = — 1- This constraint is still not suffi- 
cient, as the number of bonds at each site may still 
be equal to one or three. To get rid of the possi- 
bility of three bonds, we add a term to the Hamil- 



G 



tonian, which minimizes the number of bonds at 
each site —Jafj, where — J is the large energetic 
gain of having no dimer on a given link (J > 0). 
Therefore the full formulation of the problem is in 
terms of the Hamiltonian: 

H = -K J2 - J°fj (4) 

p(ijkl) 

with the constraint on the gauge charges Qi = 
rij=i a fj = — 1- This is equivalent to the formu- 
lation of a Z 2 gauge theory [19]. Note the, rather 
vague at this stage, analogy with the electrody- 
namics Hamiltonian H = B 2 + E 2 , supplemented 
by V • E = 0. 

Since the gauge charge operator and the Hamil- 
tonian commute [Qi,H] = 0, one can simultane- 
ously diagonalize Qi and H. However not all eigen- 
statcs H\E) = E\E) of the Hamiltonian are physi- 
cally relevant states, as they do not satisfy the odd- 
bond-per-site constraint. The projection operator 
P = Yli j(l ~ Qi) constructs a state E, that obeys 
the bond constraint, out of the cigenstates E of 
the Hamiltonian. It is also worth noting that the 
Hamiltonian is invariant under a gauge transfor- 
mation afj — ► £i<rfj€j with = ±1. Note the anal- 
ogy with the gauge transformation A — > A + V</> 
in electromagnetism. The physical observables, the 
bond field cx? (the "electric field" in the analogy) 
and the plaquette flip F p ("magnetic flux") are 
also gauge invariant. However the gauge field it- 
self is not invariant under the gauge transforma- 
tion. We ought to speak in our case of a gauge re- 
dudancy rather than gauge invariance, since the 
states related by a gauge transformation are phys- 
ically equivalent. The total Hilbert space for a N 
site square lattice model with 2N links has 2 2N 
states. Each gauge inequivalent class, however, has 
a redudancy of 2^. Therefore the number of phys- 
ical distinct states is -^tj- = 2^, which correspond 
to the fluxes through the iV plaquettes F p = ±1 . 

We can now gain insight about the confined VBS 
and deconfined RVB liquid phases of the dimer 
model from the known phase diagram of a Z 2 gauge 
theory, which is depicted in Fig. 6 in terms of the 
parameters J and K . In the J = limit, the pla- 
quette flux is fixed f p « ], as such a flux configu- 
ration minimizes the first term in the Hamiltonian 



"deconfined" (J/K) c "confined" 

Fig. 6. Phase diagram of the Zi gauge theory model Eq. (4). 
For small values of J, the plaquette flux is fixed and 
the spinons are deconfined in the spin liquid phase. For 
J S> K, the bond field is fixed, which corresponds to con- 
fined spinons in the VBS state. 

Eq.(4). The conjugate bond field afj, however, is 
strongly fluctuating. This corresponds to the res- 
onating valence bonds in the deconfined spin liquid 
state. In the opposite limit, J S> K, the bond field 
is fixed and the plaquette flux fluctuates, which de- 
scribes the behaviour of the VBS. The dimer model 
was formulated assuming the energy cost J of hav- 
ing a bond at a given link is large compared to 
K, to fulfill the one-bond-per-site constraint. From 
the phase diagram in Fig. 6, we see that for the 
quantum dimer model on a square lattice, this set 
of parameters corresponds to the confined state. In 
order for the system to be in the deconfined phase, 
one would have to relax this constraint on the value 
of J or the gauge charge constraint. The latter sce- 
nario is realized, for example, for non-bipartite lat- 
tices, such as the triangular lattice for which a RVB 
liquid phase exists [10]. 

Many things are known about the Z 2 gauge the- 
ory, which is dual to the 2D Ising model in a trans- 
verse field. In particular, we have a way to char- 
acterize the different phases through the Wilson 
loop operator Wl, defined as the product of the 
gauge operators a z around a loop of size L: Wl = 
YlieL a i- The expectation value of the operator for 
a loop of size L decays as an exponential of the 
perimeter of the loop < Wl >— e~ cL in the de- 
confined phase and as the exponential of the area 
substanded by the loop < Wl >= e~ cL in the 
confined phase. We refer to Ref. [19] for an intu- 
itive derivation of these two different behaviours. 

3.3. Analogy with electromagnetism, vison 
excitations and ground-state degeneracy 

The vaguely sketched analogy between the Z 2 
gauge theory of spin liquids and electrodynamics 
may be further exploited. It can easily be shown, 
that the bond field and the plaquette flux do not 
commute whenever the bond surrounds the pla- 
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Fig. 7. The two degenerate states in the deconfined (RVB) phase on a cylindrical topology, (a) No vison is present, (b) A 
vison has threaded the cylinder. All the links that cross the dashed line have cr? = — 1 (whereas = 1 for the rest of the 
system). 



quette. Based on the Hamiltonian Eq.(4), we can 
identify the first with the electric field and the lat- 
ter with the magnetic flux. Being non-commuting 
fields, the Heisenberg uncertainty principle applies 
and one field has to fluctuate if the other one is 
fixed. 

In the deconfined spin liquid phase, the flux 
through all the plaquettes in the ground state is 
constant (F p fa J) and the bond field fluctuates a 
lot (the dimers are disordered), as argued above. 
If the flux through one plaquette is — 1 , an excita- 
tion is created in the liquid, which has an energy 
cost of the order of K . This gapped excitation is 
called a vison and takes its name from an Ising 
vortex. Two visons on the same plaquette "anni- 
hilate" since the flux through the plaquette can 
only be ±1 - hence the reference to Ising. The 
vison is a non-magnetic (S = 0) excitation, of a 
different kind than the spinons described earlier. 
We are describing a liquid phase and both visons 
and spinons are finite energy excitations. How can 
we reconcile this with the Hastings theorem about 
the existence of low energy excitations, which tend 
to zero in the thermodynamic limit? The existence 
of the visons allows us to understand how this 
theorem may be satisfied. This theorem is fulfilled 
because the ground state is topologically degen- 
erate, when the model is defined on nonsimply 
connected spaces such as a torus. In conventional 
paramagnets, the properties such as ground-state 
degeneracy are not affected in the thermodynamic 
limit by the choice of boundary conditions. This 



is not true anymore in a topological spin liquid. 

To see this, let us take the cylinder topology, 
where a 2d system has periodic boundary condi- 
tions in one direction and free in the other. In 
the deconfined phase, locally the flux through each 
plaquette is F p w 1 , but what about the hole in 
the cylinder? Consider the flux through any curve 
C, that encircles the cylinder F^oie — ^i^co't 
(Fig. 7). In the deconfined phase of the gauge the- 
ory, we can have both situations where there is a 
flux through the hole, or no flux through it. In the 
limit K — » oo, a simple ground-state is obtained 
when all the link fields are positive of = 1. There 
is no vison present and F^oie — 1 (Fig- 7 (a)). Now 
if we flip to of = —1 all the dimers along a col- 
umn of horizontal bonds that runs the length of the 
cylinder, we have a new state where Fhoie = —1 
and where a vison has "threaded the hole of the 
cylinder" (Fig. 7 (b)). On a finite system, this is 
an excited state when J ^ 0, but in the thermo- 
dynamic limit the energy cost of creating this line 
of defect tends to zero and this state also becomes 
a ground state. Hence the two-fold degeneracy of 
the ground state. The multiplicity of the degener- 
acy of the ground state therefore depends on the 
topology of the space of system; more specifically 
on its genus in our example. For example, the de- 
generacy is four-fold on a torus. Ref. [20] describes 
in great detail the possible interpretations of this 
degeneracy and its topological nature. 

Moreover, from this example we see the hidden 
topological order is a global property of the whole 



8 



system and cannot be seen on the local scale (F p — 
1 everywhere in the "bulk" for both degenerate 
ground states) . It is easy to guess that it is impos- 
sible to deduce whether a vison has threaded or 
not the hole of the torus by looking locally at the 
dimcrs. Therefore one cannot distinguish between 
the different ground state wave functions by any 
local probe or order parameter; they are topologi- 
cally distinct. 

In the VBS phase, the bond (electric) field is 
fixed (the dimers order) and the magnetic flux 
fluctuates. This means that visons have prolifer- 
ated, they are condensed in the VBS. A site with a 
spinon has no connecting valence bonds and since 
the gauge charges are defined in terms of the num- 
ber of bonds in the system, such a spinon carries 
a gauge charge Qi = —1, which can be thought of 
as an "electric" charge. These electric charges can- 
not propagate through the fluctuating "magnetic" 
fluxes. This is another way to see that the spinons 
are confined and no longer present as finite excita- 
tions in the valence bond solid phase. 

With this electromagnetic analogy, we can also 
understand the nature of the interaction between 
the different excitations. A vison can be though 
of as a "magnetic" flux. Hence if a vison is taken 
around a spinon (the electric charge), analogously 
to the Aharonov-Bohm effect, it will result in a sign 
change in the wavefunction of the spinon ^ s — ► 
— * s . This may be thought of as a long-range "sta- 
tistical" interaction mediated by the gauge field. 

3.4. Towards realistic models 

We have shown how by melting a VBS into a 
RVB liquid one obtains deconfined finite energy 
excitations - spinons. These can be described in 
terms of a Z^ gauge theory. Exploiting the anal- 
ogy with electromagnetism, one can identify an- 
other type of excitation - visons, which correspond 
to magnetic fluxes. By considering the topology 
of the space, on which the model is defined, one 
discovers the ground states are degenerate. Hence 
Hastings' theorem is satisfied, and exotic quantum 
ground states exist. Having identified the finite en- 
ergy excitations, one can complete the electromag- 
netic analogy by considering the Aharonov-Bohm 



type interactions between the "electric" and "mag- 
netic" charges. 

Here we have only focused on phenomenological 
description of spin systems which can have a topo- 
logical Zi spin liquid phase. Actually, there are 
already several models which can be shown (some- 
times on rigorous grounds) to sustain Zi topo- 
logical order. Without being exhaustive, one can 
mention quantum dimcr models on non-bipartite 
lattices [10], Ising-like models with multiple spin 
interactions [11], rotor bosonic models [12] and 
more recently models Heisenberg-like spin models 
with original SU (2) symmetry [13] or without [14]. 
We should also note that the Z^ spin liquid is 
only the tip of the iceberg: there are many much 
more intricate topologically ordered phases pos- 
sible, some with fractional and even non-Abclian 
statistics [15]. These phases are very promising 
candidates as possible realizations of quantum 
computers, as they could be topologically pro- 
tected from decoherence effets [16]. We also have 
kept the discussion to 2d systems, but this is not 
essential to obtain topological order. It is, how- 
ever, necessary to obtain the fractional statistics 
and non-abelian properties. 

Finally, we note that at this point theoretical 
models are sadly not very realistic. It remains a 
huge challenge to develop theoretical techniques to 
look for topological spin liquids in realistic mod- 
els and find them in the laboratory. There is no 
experimental evidence, except for the fractional 
quantum Hall effect, as to the existence of topo- 
logically ordered phases. It was theoretically pro- 
posed [17] that topological order and visons could 
exist in high T c cuprates superconductors. Later 
experiments [18] disproved this statement. Never- 
theless, these ideas and the proposed experimental 
set-up to detect visons could be useful to eventu- 
ally experimentally engineer a real topological or- 
dered phase, for example in a Josephson Junction 
array. The key ingredient is certainly to try to dis- 
order the system as much as possible. 
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4. Algebraic or Critical Spin Liquids 

4.1. Description of algebraic spin liquids 

Besides topological degeneracy, an alternate way 
in which an apparently featureless quantum para- 
magnet can remain "gapless" as required by the 
Hastings' theorem, is by having gapless bulk ex- 
citations - i.e. being a "critical" quantum phase. 
Such critical quantum ground states have power- 
law forms for correlation functions, but no true 
long-range order. The lack of a finite order param- 
eter is necessary for a phase to be called a spin 
liquid. In some sense, these systems behave as if 
they were at a critical point without being tuned. 
A physical picture in terms of the spin degrees of 
freedom for such critical spin liquid phases can be 
visualized in terms of valence bonds between pairs 
of spins separated by arbitrary distances. The va- 
lence bond lengths follow a power-law distribution. 
These valence bonds quantum-mechanically fluc- 
tuate, which produces the continuum of excitations 
mentioned above. The spinons (unpaired spins), 
unlike in topological spin-liquids, arc strongly in- 
teracting with this valence bond background and 
cannot be considered as free particles. A cartoon 
picture of a critical spin liquid is given in Fig. 3c. 

Much less is known about critical spin liquids 
than topological spin liquids, as introduced in 
Sec. 3. Technically, the theoretical models in- 
volved are more complicated and even toy models 
are difficult to construct. Apart from the general 
considerations given in the preceding paragraph, 
there is no clear overall theoretical picture of all 
the generic features of critical spin liquids, al- 
though some progress has recently been made in 
this direction [22]. However, an advantage over 
topological spin liquids is that there are several 
experimental spin systems which are good candi- 
dates for being a critical spin liquid, and can thus 
be confronted with theory. We will describe one 
particular theoretical scenario which gives rise to 
a critical spin-liquid phase, which is consistent 
with current experiments on one promising com- 
pound. This example will let us identify some of 
the elements that are needed to form a critical 
spin-liquid. 



4.2. Triangular lattice S = 1/2 antiferromagnets 
- Experimental evidence of spin liquid behaviour 

As already mentioned, the key idea to obtain a 
spin liquid is to enhance the conditions that can su- 
press magnetic order, namely using low spins (S = 
1/2 or S = 1), low dimensional systems, and geo- 
metric frustration, as in non bipartite-lattices [23] . 
The Id Heisenberg chain is the simplest example 
of a critical system - it is much harder to find one 
in 2d. Our interest lies in 2d systems, thus we will 
consider a spin 1/2 antiferromagnetic sytems on 
the triangular lattice. Numerical calculations have 
determined that the ground state of the isotropic 
Heisenberg model on the triangular lattice has long 
range coplanar magnetic order, with spins order- 
ing in a 120 degrees fashion. This order is the 
same as for the classical ground state of the model, 
which results from the minimization of the frus- 
trated interactions. However, in the S — 1/2 quan- 
tum system, the magnetization is significantly re- 
duced compared to the S — > oo classical limit, due 
to quantum fluctuations. This numerical result, 
along with recent experiments, gives one reason to 
believe that a spin liquid state is perhaps accessi- 
ble in such triangular lattice antiferromagnets. 

Currently, two experimental systems are likely 
candidates for spin liquid behaviour in the tri- 
angular S = 1/2 antifcrromagnet: an organic 
material k-(ET) 2 Cu 2 (CN) 3 (Ref. [24]) and a tran- 
sition metal material CS2CUCI4 (Ref. [25]). We 
will focus the discusion on the latter system for 
which the Cu atoms have nine electrons on the d 
shells and the one missing electron results in the 
5 = 1/2 moment. The experimental phase dia- 
gram for CS2CUCI4 (Fig. 8), shows a long range 
magnetically ordered phase at very low tempera- 
tures and a "spin liquid" in an intervening tem- 
perature regime. At higher temperatures the spins 
on the Cu atoms fluctuate independently and 
the magnetic susceptibility obeys the traditional 
Curie- Weiss law. The planar triangular lattices 
in the Cs 2 CuCLt crystal are stacked on top of 
one another, with a large interplanar distance. 
Therefore the in- plane interactions J and J', as 
depicted in Fig. 9, are dominant. The magnetic 
long range ordered state is coplanar as predicted 
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Fig. 8. Experimental results obtained on CS2CUCI4. Left: Low field magnetization versus temperature (taken from Ref. [25] d). 
At the transition from the long-range ordered (LRO) state at low temperatures to the short-range ordered (SRO) spin liquid 
at intermediate temperatures the magnetization exhibits a cusp. The broad peak at the transition from the spin liquid to 
the paramagnet present at high temperatures is characteristic of short-range antiferromagnetic correlations. Right: Neutron 
scattering intensity as a function of energy transfer (taken from Ref. [25] b). The sharp peak at low energy is the signature 
of long-range ordering. The broad continuum at high energies can be fitted by a power-law (see inset) and is ascribed to 
the low-energy excitations of a critical spin liquid state. 




Fig. 9. The long-range coplanar ordered state adopted by 
the S = 1/2 antiferromagnet Cs 2 CuCl 4 (see Ref. [25]). The 
bond directions <5i and 82 and the interaction strengths J 
and J' in the effective Hamiltonians Eq. (6) and Eq. (7) are 
also shown. 

by numerical studies. Upon moving through the 
lattice the spins spiral in the plane. However, 
since the exchange interactions J and J' differ in 
magnitude the wavelength of the spiraling spins is 
incommensurate with the lattice, resulting in the 
spin ordering represented in Fig. 9. 

In CS2CUCI4 the experimentally detected "spin 
liquid" only exists at finite temperature, whereas a 
theoretical spin liquid is well defined only at T = 0. 



Despite this drawback, there are two significant ad- 
vantages of CS2CUCI4 that compensate. Firstly, it 
is possible to grow large high purity single crystals, 
which may be used in neutron scattering. The main 
quantity which can be accessed in neutron scatter- 
ing experiments is the magnetic structure factor, 
which is the Fourier transform of the spin-spin cor- 
relation function. At low temperatures, the exper- 
iments show magnetic Bragg peaks at wavevectors 
corresponding to the wavelength of magnetic or- 
dering, which detect the spiral long-range ordered 
state (see Fig. 8). Inelastic neutron scattering ex- 
periments at higher temperatures, where the long 
range order is destroyed, can measure the strength 
of quantum fluctuations in the putative "spin liq- 
uid" . Strikingly, these measurements show a power 
law dependence on energy transfer u) of the struc- 
ture factor 

(S + (q,cj)S-(-q,-u)) ~ c^ 2 (5) 

with 77 ~ 0.75 for wavevectors q chosen close to 
those of the Bragg peaks. This power-law feature 
is suggestive of a critical spin-liquid phase. The in- 
tensity of the scattered beam as a function of lu is 
plotted in Fig. 8b. Two features can be observed. 
Firstly, a sharp low-energy peak which shows up 
at low temperature is the signal of spin-waves as- 
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sociated with the long-range order. This peak dis- 
appears for T > T c . Secondly, a broad continuum 
scattering which extends up to high to suggests 
many low-energy excitations. The number of these 
low-energy excitations increases as the aforemen- 
tioned power law, as the transfer energy is lowered 
(see log- log scale in the inset of the figure) . 

The second advantage of CS2CUCI4 is that the 
exchange interactions J and J' are very small, of 
the order of a few Kelvins. Therefore one can exper- 
imentally measure the parameters of the effective 
Hamiltonian. These experiments are conducted in 
high magnetic fields, where the spins are forced to 
be polarized along the external field. This is possi- 
ble because the Zeeman energies involved are of the 
order of a few Kelvins, corresponding to magnetic 
fields of a few Teslas, which can easily be reached 
in the laboratory. Using a neutron beam a single 
spin may be flipped and by measuring the energy- 
momentum dispersion relation for this flipped spin, 
the values of the exchange interaction may be ex- 
tracted. The Hamiltonian for a single layer has the 
form of nearest neighbour exchange interactions on 
the triangular lattice (see Fig. 9) 

Ho = ^[JS r • S r+ 5 1+ 5 2 

r 

+ J'(S r • Sr+Sj + S r • S r+ 5 2 )]. (6) 

In addition there is a term due to spin-orbit 
coupling, called the Dzyaloshinskii-Moriya inter- 
action. In most systems this term is small, but it 
may not be neglected here, since the exchange in- 
teractions are small. The full Hamiltonian there- 
fore has the form H = H + Hdm, where the 
Dzyaloshinskii-Moriya interaction takes the form: 

Hdm = - ^2 D • S r x (S^^ + S r+< 5 2 ) (7) 

r 

where the experimentally measured parameters 
are J - 4.3K, J'/J - 0.34 and \D\/J ~ 0.053 
(Ref. [25]c). The Dzyaloshinskii-Moriya interac- 
tion breaks the SU(2) symmetry down to a U(l) 
symmetry, called the "easy-plane" symmetry. 

The model presented cannot be solved analyti- 
cally or numerically except on very small lattices. 
Again due to the frustrated interactions, Monte 
Carlo calculations are plagued by the sign prob- 



lem. To gain an understanding of the physics in 
the particular phases, it is instructive to consider 
a simpler model, which has the same symmetries 
as the Hamiltonian in Eq. (6) and (7). The fact 
that the Dzyaloshinskii-Moriya interaction breaks 
the SU{2) to U{1) symmetry suggests to use an 
anisotropic exchange interaction Jf- < Jij , leading 
to the following easy-plane spin model: 

Hxy = YM^ S 1 + S ? S l) + (8) 

where the notation reminds us that there are 
two different coupling strenghts J and J' in the 
real material CS2CUCI4. 

4.3. Chain of transformations: from spins to 
bosons, from bosons to bosonic vortices and from 
bosonic to fermionic vortices 

Even this simpler model is still hard to solve. 
The theoretical challenge is to find other vari- 
ables, which correctly describe these spin systems. 
Here we present one of the proposed approaches 
[26], which allows one to make progress, based on 
various transformations. In the limit J?- = 0, the 
model in Eq. (8) is the quantum version of the 
XY model. In 2d, we know that the classical XY 
model can be equivalently described in terms of 
vortices. An analogous duality transformation in 
the quantum system will lead to the description 
of the system in terms of quantum vortices. A 
second transformation, which fermionizes the vor- 
tices, which are bosons, proves useful to make the 
system tractable. This change of variables, called 
the Chern-Simons transformation [27], is the 2d 
generalization of the well-known Jordan- Wigncr 
transform in Id. 

First, we conveniently describe spin 1/2 objects 
in terms of hard-core bosons. A spin up on a given 
site corresponds to the site being occupied by a 
boson, whereas a spin down corresponds to an un- 
occupied site. Hard-core bosons means that there 
is at most one boson per site. Then the x — x 
and y — y spin interactions are mapped onto the 
hopping of bosons, which is described by creation 
Sf = Sf + iSf and annihilation operators 5~ = 
Sf — iSf . On a bipartite lattice, the sign of spin in- 



12 



teractions (ferromagnetic or antiferromagnetic) is 
irrelevant, since the sign of the hopping term can 
be easily removed with a simple unitary transform. 
But for non-bipartite lattices, such as the triangu- 
lar lattice, antiferromagnetic interactions between 
the spins lead to a positive sign for the hopping 
term. Within a boson picture this does not make 
sense since the bosons energy must be lowered by a 
near-neighbor hop. In order for the kinetic energy 
to have the "correct" negative sign for the hopping 
strength, one needs to do a canonical transforma- 
tion to include a fictitious flux of tt through the tri- 
angles of the lattice. When the boson hops around 
a triangular plaquette, its phase picks up a factor 
of 7r, reestablishing the sign due to antiferromag- 
netic interactions. This tt flux produces a fictitious 
magnetic field, in which the boson hopping takes 
place. The Hamiltonian Eq.(8) in terms of bosons 
then takes the following form 

H *v = -Y,[^f S tS;e iA % + h.c. + Jf^S]}, (9) 

m 

where the last term is interpreted as a nearest 
neighbour density-density interaction between the 
bosons and where [, . , A • dl = tt. 

J triangle 

We now proceed with the next transformation. 
Bosons (in magnetic fields) are often described in 
terms of vortices [28] . Let us now focus on this new 
object and its properties. To define a vortex, one 
can consider the phase of the boson creation oper- 
ator ~ e 1 ^ (or equivalently the spin direction 
in the XY plane). A vortex is present when the in- 
tegral of the gradient of the phase <fi around a loop 
gives J V</> • dl = 2ir. Depending on the configura- 
tions of the spins in the triangle, they may result in 
a vortex, associated with a vortex number N = 1, 
or an antivortex, associated with a vortex number 
N = 0, as shown in Fig. 10a. 

Since the original spin problem is frustrated, the 
two spin configurations in Fig. 10a are energetically 
equivalent. Therefore there are as many triangular 
plaquettes with vortices as with antivortices; the 
dual vortices are at half-filling. 

In this description, the vortices live on the hon- 
eycomb lattice, which is dual to the triangular lat- 
tice. This is easy to see, by connecting sites at 
the centers of plaquettes of the triangular lattice - 




Fig. 10. (a) A pictorial representation of vortices and an- 
tivortices. (b) The honeycomb lattice is dual to the trian- 
gular lattice. 

see Fig. 10b. The duality transformation between 
bosons (the original spins 1/2) and vortices intro- 
duces a gauge field a sitting on the links of the 
dual lattice. This gauge field will encode the long- 
range interaction between the vortices. This exact 
mapping relates, for example, the z component of 
the spin to the field: S z + 1/2 = x a|. This 

indicates that a vortex hopping around a dual pla- 
quette will pick up a tt phase factor. This is easy 
to understand as a vortex hopping around a pla- 
quette of the honeycomb lattice will encircle one 
spin 1/2, and thus "feel" its average flux. Since on 
average (S z ) = (the original spins are not in a 
magnetic field) , the vortices experience an average 
background magnetic field of flux J d 2 rV x a = 
tt. a is the static part of the gauge field, precisely 
designed to account for the fact that the original 
spins can take only the values S z = ±1/2. 

The gauge field actually mediates a long-range 
Coulomb, or "electromagnetic", interaction be- 
tween the vortices. In 2d, the Coulomb potential 
is logarithmic. The Hamiltonian in the vortex rep- 
resentation, where and b are the vortex creation 
and annihilation operators, reads: 

H = - ^{Ujblbje^+O + h.c.) 

(ij) 

+ Y,{N i -l/2)V ij {N j -l/2) 

ij 

+ [/]T(Vxa)2. (10) 

r 

The first term is the kinetic energy of the vor- 
tices, moving on a honeycomb lattice in a field 



13 



with static component satisfying (V x a°) r = n 
for all sites r of the triangular lattice. The second 
term is the Coulomb interaction between vortices 
with Vij ~ ln(|i — In this term we clearly see 
the half-filled natures of the vortices, as Ni = b\bi 
is the vortex number operator. The final term - 
similar to a magnetic energy term - accounts for 
the dynamics of the gauge field. We refer again to 
Rcf. [26,28] for a rigorous construction of the dual 
vortex hamiltonian. 

This type of boson/vortex transformation is es- 
pecially useful when the dual vortex variables are 
at integer filling. Here, however, due to the original 
spin frustration, the vortices are located at half- 
filling, since the states with filling fraction of zero 
or one are equally likely. The dual vortex repre- 
sentation thus appears as intractable as the origi- 
nal spin model. However, further progress can be 
made by performing yet another transformation, 
"transmuting" the bosonic vortices into fermionic 
vortices. The main advantage of this change of 
variables, is that the interacting boson field may 
be described using a non-interacting fermion field, 
since fermion statistics have Pauli exclusion built 
in. This makes the model mathematically easier. 
However, there is a price to pay. The fermion wave- 
function changes sign upon exchanging particles, 
whereas the boson function does not. The transfor- 
mation must therefore correct for this sign change. 
This Chcrn-Simons transformation [27] may be in- 
tuitively understood, as attaching a fictitious 27r 
flux tube to a fermion vortex in order to represent 
the bosonic vortex. So if we consider the exchange 
interaction between two fermionic vortices, due to 
the Aharonov-Bohm effect, encircling one particle 
half way around the other results in a e lTr flux from 
the interaction of the fermion with the flux tube of 
the other particle, which cancels the negative sign 
from fermion statistics. 

Now let's look at a static picture of the fluxes 
attached to these fermions, which also live on the 
dual lattice. We can use a convention where each 
fermion shares on average its 27r flux amongst the 
three honeycomb plaquettes to which it belongs. 
Since on average, each honeycomb plaquette has 
three vortices, as the vortices are at half-filling, it 
is therefore pierced by a 2n flux. This is equivalent 



to a zero flux and therefore irrelevant for the be- 
haviour of the fermions. Of course, this is a crude 
picture, but it is valid in the mean field approxima- 
tion where we ignore flux fluctuations around their 
average. Therefore, the fluxes associated with the 
Cherns-Simon transformation cancel to zero when 
we smear out these fluxes in a mean field fashion. 
We note that any mean-field approximation is of 
course crude, but here it is actually expected to 
be rather good, since the interactions between the 
vortices tend to supress their density fluctuations, 
and therefore the flux fluctuations in the Cherns- 
Simon language. 

So now in the mean-field picture we are left 
with free fermions. However, we have to remem- 
ber that these fermions are also coupled to the 
original gauge field a, as were the vortices. In 
mean-field theory, we only consider the static con- 
tribution from a , and we are thus describing a 
system of free fermions on the honeycomb lattice 
with a background n flux. This is a one-particle 
problem and it is easily solved. The corresponding 
Fermi "surface" consists of only four points, with 
Dirac-like (gapless) dispersions around each one. 

The non-interacting limit we are currently 
describing, is a good starting point to now re- 
introduce the gauge fluctuations we supressed in 
the mean-field approximation. To do this, we now 
describe the behaviour of the low energy particle- 
hole excitations around these four Fermi points 
in terms of a Dirac theory for fermionized vor- 
tices. There are four flavours of Dirac particles, 
which correspond to the four Fermi points. These 
fermionic vortices interact logarithmically, which 
is accounted for by a coupling to the gauge field 
a. Additionally, they are coupled to the Cherns- 
Simon gauge field Acs due to the flux tube. It 
can be argued [26] that the field due to the flux 
tubes fluctuates much more than the field due to 
the Coulomb interactions between vortices and 
that at long length scales the Chern-Simons field 
can be ignored. This may also be seen explicitly 
by defining an auxilary field and integrating out 
the field due to the fictitious flux tubes. In fact, 
the logarithmic interaction in 2d is so strong that 
it dominates over the statistics encoded by the 
Cherns-Simon flux tubes. In this respect the de- 
scription in terms of fermions and bosons is anal- 
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ogous. We invite the interested reader to consult 
Ref. [26] for a more detailed discussion. 

The final formulation of the critical spin liquid 
in terms of N — 4 flavours of Dirac particles, which 
interact by a logarithmic potential in 2d is in terms 
of a QED (Quantum Electrodynamics) Lagrangian 
in 2 + 1 dimensions 

Lqeds = ^«7"(fy ~ + ^ (V x a) 2 , (11) 

where ip a is the spinor associated to the Dirac 
fermion flavour a e {1,2,3,4}, are Dirac 
gamma matrices in direction /i G {x,y,z} and 
a = a + Acs is the auxiliary field used in the inte- 
gration of the Cherns-Simon field ( the sums over 
double indices are implicit). Strictly speaking, 
there should be an additional L^f four- fermion 
term that originates from short-range interactions 
between the vortices that we have not taken into 
account explicitly, but that can always be present 
in the microscopic model. It is known that such 
terms are irrelevant in the limit of an infinite 
fermion flavour number N — ► oo and for all N > 
N c , where the exact value of N c is not known. 

Assuming N c < 4, we can safely ignore this 
term. Now, we can finally import the results from 
QED, which indicate that Eq. 11 describes a crit- 
ical theory with gapless excitations. Thus as a re- 
sult of these consecutive changes of variables, from 
spins to Dirac particles, this approach describes 
the much desired critical spin liquid phase, with- 
out fine tuning of variables and no free particles. 
This specific spin liquid is dubbed an "algebraic 
vortex liquid" [26] . 

In this section, we have tried to describe the 
different steps of the construction that lead to the 
algebraic vortex liquid. We will now only men- 
tion some of the properties of this type of spin 
liquid. The critical phase not only does not break 
any symmetries, but a 5J7(4) symmetry emerges. 
Above all, the theory is able to make predictions, 
which may be compared with experiment. The 
structure factor is predicted to have power law 
correlations around five distinct momenta with 
the same exponent for each wavevector. Two of 
these momenta correspond to the experimentally 
measured magnetic long range order wave vectors. 
To make a further comparison with CS2CUCI4, 



the Dzyaloshinskii-Moriya interaction be incor- 
porated, and in the Dirac theory corresponds to 
a "mass" term which destabilizes the algebraic 
vortex liquid driving it into an incommensurate 
spiral state at low enough temperature: this is ex- 
actly what is observed in the experimental phase 
diagram of CS2CUCI4 (see Fig. 8a). 

The proposed algebraic vortex liquid state is not 
the only existing scenario to describe the physics 
of CS2CUCI4. There have been recent theoretical 
proposals such as a quasi-ld approach [29], an alge- 
braic spin liquid state with fermionic spinons and 
a U(l) compact gauge field [30] or the proximity of 
a quantum critical point between a Z 2 spin liquid 
and a spiral phase [31]. 

We would finally like to emphasize that through 
the example described in this section, we can see 
that, like in the topological spin liquid case, the 
natural description of these phases is in terms of 
gauge fields. The gauge fields that emerge for criti- 
cal spin liquids, however, have different symmetries 
(here C/(l)) than the Zi symmetry for the topolog- 
ical spin liquids of Sec. 3. We will further exploit 
gauge theories in the next section. 



5. Deconfined quantum critical points 

We will now turn from the description of ex- 
otic phases of matter discussed in the two previ- 
ous sections, to an exotic critical point separating 
two conventional states of matter. We will focus 
on the quantum phase transition between a Neel 
state, which breaks spin rotational symmetry, and 
a Valence Bond Solid state, which breaks trans- 
lational symmetry, depicted in Fig. 1 and 2b. We 
will show that, contrary to general expectations of 
the Landau theory of phase transitions, a direct 
second order transition is possible. It is a contin- 
uous phase transition between two ordered states 
with different symmetries. We will find that, right 
at the critical point, the spinons we encountered in 
the two previous sections are essentially deconfined, 
whereas they are altogether absent in both sur- 
rounding phases. The description of this so-called 
"deconfined quantum critical point", will also be 
in terms of a gauge theory, namely a non-compact 
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Fig. 11. The T = transition (driven by an external parameter g) between a Neel antiferromagnetic order and a Valence 
Bond Solid for spins S = 1/2 on the square lattice. The transition is argued to be second-order [32] with both Neel and 
VBS order parameters vanishing continuously at the critical point. 



CP 1 theory in this case. 

5.1. Formulation of the problem: the transition 
and standard Landau theory 

We want to describe 2c? quantum S =1/2 spin 
systems with one spin per unit cell on the square 
lattice and the quantum phase transition, i.e. a 
transition at zero temperature driven by an exter- 
nal parameter, between a Neel state and a VBS 
state (see Fig. 11). We do not explicitely men- 
tion a specific form of the microscopic interaction 
strength g that drives the transition, which could 
be longer-range couplings like in Hamiltonian 
Eq. (3), or maybe a ring-exchange term, but rather 
focus on the symmetries of the associated phases. 

If instead of a VBS state a simple paramagnet 
with no symmetry breaking existed on the right- 
hand side of Fig. 11, it would be natural to de- 
scribe the transition within a standard Landau 
approach. That is, we would identify a local or- 
der parameter, the staggered magnetization n,; = 
(— l) ri Si, coarse-grain it to the continuum n(r) and 
expand the free energy or Lagrangian in gradients 
or powers of n(r) allowed by the symmetries such 
as i(n) = |Vn| 2 + r|n| 2 + u|n| 4 + . . .. Depending 
on the precise value of the parameters, we would 
get either the Neel state (r < 0) or the paramagnet 
(r > 0), and the transition at r = 0. Here the main 
problem is that the plain paramagnet does not ex- 
ist in 2D at T = 0, according to Hastings' theorem. 
It is thus natural to consider something else than 
the paramagnet, more specifically the VBS. In the 
case of a transition between two different broken 



symmetry phases, the Landau theory would ex- 
pand the Lagrangian in terms of both the Neel or- 
der parameter n(r) and the VBS order parameter 
V'vbs with possible cross terms. 

For a phase transition between states with dif- 
ferent broken symmetries, the resulting Landau 
theory predicts one of the four possible scenarios: 
(i) a direct first order phase transition, (ii) an in- 
termediate coexistence phase, between Neel and 
VBS orders in our case, (iii) an intermediate phase 
with no order at all, (iv) a second order transi- 
tion in the case of a multi-critical point. Scenario 
(iii) can be excluded using Hastings' theorem. Sce- 
nario (iv) is very unlikely, because it is not generic 
and it implies a fine-tuning of some parameter. 
Scenarios (i) and (ii) cannot be ruled out on gen- 
eral grounds, but both exclude a direct continuous 
(second-order) transition. 

In contradiction to this standard Landau analy- 
sis, recent work [32] shows that in the specific case 
of a Neel- VBS transition on the square lattice, the 
transition can generically be second order without 
any fine-tuning, due to subtle quantum interfer- 
ence effets that invalide the Landau analysis. The 
crucial elements responsible for a continuous tran- 
sition are topological objects called hedgehogs. The 
resulting field theory is a non-compact gauge the- 
ory. In what follows, we review in a simplified way 
the approach of Ref . [32] . 

5.2. Neel order, skyrmions and hedgehogs 

To gain insight into such a Neel- VBS transition, 
let us first consider the Neel state, and its topologi- 
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Fig. 12. Real-space representation of a skyrmion in the Neel field n at a given timeslice. Spins in the center are pointing 
down whereas spins on the boundaries point up. The charge of this skyrmion is Q = +1. 



cal defects. Since we are considering a 2d quantum 
system at T = 0, the Neel vector lives in (2 + 1) 
dimensions, where the extra-dimension can be in- 
terpreted as imaginary time in the standard path 
integral formalism. We shall denote the space-time 
coordinates as x^ with fi 6 {x,y,r}. Consider 
smooth configurations of the staggered magnetiza- 
tion n(Xft) = (—l) x+v S(x fl ), i.e. n(x^) is slowly 
varying and is of constant amplitude. Let's define 
a unit length Neel vector n(x M ) = n(a; AI )/|n(x A1 )|. 
For the Neel ground states, all the n(x M ) are point- 
ing in the same direction, but there exist configu- 
rations which admit a topological defect called a 
skyrmion. A skyrmion is represented in Fig. 12, in 
which the Neel vector at the center of the "sample" 
points down whereas it points up at the borders. 
Mathematically, the skyrmion number at a given 
time slice can be defined as an integer topological 
number 

Q = — I dxdy(n ■ d x ii x cLn) (12) 

47T J 

sometimes called an integer charge. If the spins 
were living in the continuum, this number would be 
a conserved quantity in time. However, on a lattice, 
this number can change from one time slice to the 
other. Such a space-time event where the skyrmion 
number Q changes by ±1 is called a hedgehog (see 
Fig. 13a). The choice of this name hopefully be- 



comes apparent upon inspecting Fig. 13b, where 
a hedgehog is represented. We can clearly see a 
singular point from which all the n are pointing 
outwards, like quills on a hedgehog. A hedgehog is 
the term used for a configuration of unit Neel vec- 
tors, which is singular at one space-time point, but 
smooth elsewhere. What is the role of these topo- 
logical defects? In a Neel state, hedgehogs are very 
costly energetically and are therefore absent. Deep 
within an ordinary paramagnet, the spins fluctu- 
ate essentially independently of one another. In 
this case there is no energetic reason why hedge- 
hog events should be suppressed: we say that the 
hedgehogs proliferate. 

To see what happens to these hedgehogs on the 
VBS side and at the critical point, it is useful to 
consider the non-linear sigma model description of 
the Neel state [33] with the following action: 

Sa = j^Jdx ll \d ll n\ 2 + S Beiiy (13) 

where the first term accounts for the slow varia- 
tions of n, g is some "stiffness" constant, and the 
second term is a Berry phase term which accounts 
for the quantum nature of the S = 1/2 spins. Now 
let's see the effects of hedgehogs on such a formu- 
lation by expanding the partition function Z = 
J dnexp(-Sa) in powers of the fugacity A of such 
events: 
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Fig. 13. (a) The skyrmionic charge is suddenly changed at a hedgehog event in space time, (b) Real-space representation 
of a hedgehog event (the two spin configurations represent different time slices). A hedgehog corresponds to a singular 
configuration of n(x^) at one space-time point where the skyrmionic charge changes. All spins are pointing outwards of 
the hedgehog. 



z = z + j KxDz^xl) 

'I 

+ \ J X(xl)X(xl)Z 2 (xl,xl,) + ...(U) 

where Zq describes a model where no hedgehogs 
are allowed. It can be shown (see Ref. [33] a) that 
due to the Berry phase terms on the square lat- 
tice, the contribution Zi(xV) from single hedgehog 
event is oscillatory in space-time and thus rapidly 
vanishing upon integration over xz - it is there- 
fore strongly "irrelevant" . The same is true for Z2 
and Z3, and the first non- vanishing contribution 
is from quadruple-hedgehog events (A4 ~ A 4 ) cor- 
responding to skyrmion number changes of ±4. In 
powers of the one-hedgehog fugacity, this is already 
a "rare" event, as it requires four hedgehog events 
acting at the space-time point. It is compellingly 
argued in Ref. [32] that A4 is also irrelevant in a 
renormalization group (RG) analysis right at the 
critical point. Therefore, hedgehogs are completely 
absent, in a RG sense, at the phase transition. 

Now what is the role of the hedgehogs in the 
VBS? Whereas they are completely irrelevant in 
the Neel state, hedgehogs turn out to be crucial 
for the onset of VBS order. It can be rigorously 



shown [34] that the skyrmion creation operator, 
identified with the presence of a hedgehog, is equal 
to the VBS order parameter ipVBS up to constant 
(see Ref.[32]b for an intuitive derivation). This es- 
sentially means that when hedgehogs proliferate, 
they do so in such a way that the spins pair into sin- 
glets and form a valence bond solid. This is again 
due to Berry phase effects. This dual property of 
hedgehogs - simultaneously destroying Neel and 
creating VBS orders - cannot be captured within a 
standard Landau analysis. The emergence of VBS 
order is moreover in agreement with Hastings' the- 
orem which precludes a conventional gapped para- 
magnet. 

So, the non-linear sigma model, supplemented 
with an analysis of the role of hedgehogs, contains 
in a way all the physics of both Neel and VBS 
phases. However, it is not possible to build a Lan- 
dau theory out of it to describe the transition be- 
tween both states, because hedgehogs are absent 
right at the critical point as we argued above. 

5.3. Gauge theory description 

We now consider an alternate formulation which 
enables us to glean more insight into the physics of 
the deconfined quantum critical point. Once again, 
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gauge theory will provide the natural framework. 
We first rewrite the Neel vector in terms of spinor 
fields z a , where a =|, I, 

n = z^a^pzp, (15) 

where a is the vector of Pauli matrices and the sum 
over double indices is implicit. Also, the normal- 
ization |n| 2 = 1 requires |zj 2 + |z|| 2 = 1. This is 
the so-called CP 1 representation of the SU (2) al- 
gebra. The spinors z a can roughly be interpreted 
as the spinon fields. Eq. 15 simply means that we 
need two of these S = 1/2 objects to get the S = 1 
spin-flip Neel field. 

By definition, the spinor fields have a local U(l) 
gauge redundancy - if we multiply z by some 
phase factor exp(i^(x fl )), n stays unchanged. The 
spinons are therefore coupled to a U(l) gauge field 
= ^s(z^dfiz). The non-linear sigma model ac- 
tion without Berry term can therefore equivalcntly 
be written as 

So- = ^ J dx^d^-ia^z^ 2 . (16) 

Where are the hedgehogs in this formulation? 
Within the CP 1 representation, the skyrmion 
number Q can be shown to be simply related to 
the flux of the gauge field 

Q = J dxdy(d x a y - d y a x ). (17) 

A skyrmion is interpreted as a flux charge, and 
the hedgehog which changes this number is there- 
fore a monopole in the gauge flux. A monopole 
event is sometimes called an instanton. Precisely at 
the critical point, the hedgehogs/monopoles have 
been argued to be irrelevant and their absence 
gives rise to a supplementary conservation law: 
the gauge flux is a conserved quantity. In absence 
of monopoles, a gauge theory is said to be non- 
compact. Thus a non-compact CP 1 gauge theory 
is a likely candidate to describe the physics of the 
critical point, since it incorporates from the very 
beginning both the spin 1/2 and the absence of 
hedgehogs. A Landau theory in terms of both order 
parameters of the Neel and VBS phases does not in- 
corporate this crucial feature, and thus also misses 
the important fact that the spinons, although ab- 
sent in both phases, are the good "emerging" de- 



grees of freedom to describe the system at the crit- 
ical point. 

5.4. Non-compact CP 1 theory and phase diagram 

With these arguments in hand, the authors of 
Ref. [32] now argue that the Neel- VBS quantum 
critical point can be described by a non-compact 
CP 1 gauge theory with the following Lagrangian 

^nccpi = \(du.- ia^Za] 2 + r\z\ 2 + u(\z\ 2 ) 2 

+K(e liI/ xd v ax) 2 , (18) 

where is the antisymmetric tensor, and we 
recognize a magnetic-like energy in the last term. 
This is a Landau-like expansion, but in term of 
the spinon fields, which can admit a second-order 
phase transition, at least within a mean-field anal- 
ysis. Specifically, for r < 0, the model describes the 
Neel state, since the spinon field acquires a non- 
zero expectation value (z a ) ^ and thus (n) = 
(4) a <*,P i z ») ¥= 0. For r > 0, (z a ) = (n) = 0, and 
we have a paramagnetic state, also with no hedge- 
hogs. This hedgehog-free state is not an ordinary 
gapped quantum paramagnet. Rather, it turns out 
to be a critical U(l) spin liquid, as described in 
Sec. 4. For r = 0, we have a second-order transi- 
tion between these two states: therefore the Neel- 
VBS critical point is governed by the physics of a 
Neel-f (1) spin liquid transition. 

How can we gain intuition about this new phase 
transition? Refs. [32,35] argue that it can be de- 
scribed by a simpler classical model: a 3d Heisen- 
berg 0(3) model with no hedgehog topological de- 
fects. Such a model is hard to define analytically; 
it corresponds to the notation Z in Eq. (14). How- 
ever on a computer, a Monte Carlo calculation 
can be performed by simulating a classical 0(3) 
magnet and rejecting "by hand" all configurations 
which have a hedgehog. Such a calculation was per- 
formed [35] and a second order phase transition 
was found, between the ferromagnetic state and 
the no-hedgehog paramagnet, with a different set 
of critical exponents than if one allows free hedge- 
hogs to exist. These calculations showed that a de- 
scription of the transition can be given in terms of 
the spinon fields as in Eq. (18), and not in terms of 
the Neel field. The observed transition is argued to 
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Fig. 14. Renormalization-group flow phase diagram with 
parameters A4 (the quadruple-hedgehog fugacity) and g 
(the coupling strength that drives the transition). Four 
fixed points are visible (see text): the Neel (denoted AF), 
the VBS, the [7(1) spin liquid ([7(1) SL) and the deconfined 
quantum critical point (DQCP). 



be in the same universality class as the Neel- VBS 
2d quantum phase transition. 

The "RG flow" phase diagram in terms of the 
fugacity of the quadruple-hedgehog event accounts 
for all the physical phases involved in the above 
discussion (see Fig. 14). The flow admits four fixed 
points. The first point is the antiferromagnetic Neel 
phase, where the fugacity renormalizes to zero, 
since the hedgehogs are too costly, and the VBS 
phase where hedgehogs proliferate, since they cre- 
ate the VBS order. Then a novel deconfined quan- 
tum critical point is present, where 4- hedgehogs 
events have been argued to be irrelevant. Finally, 
the phase diagram admits a U(l) spin liquid point, 
which can be argued on general grounds [36] to 
be unstable with respect to the insertion of hedge- 
hogs events. The deconfined quantum critical point 
might be hard to observe in numerical calcula- 
tions, because in a lattice model the "bare" A4 is 
in general non-zero. However, there are many ex- 
perimental predictions that can be made based on 
such a theory, in order to be able to detect decon- 
fined criticality in the laboratory. We refer the in- 
terested reader to references [32] to find these pre- 
dictions and to have a complete overview of decon- 
fined quantum criticality. Other pedagogical intro- 
ductions to this topic are also available [37] . 



6. Conclusion 

We have discussed how "exotic" behaviour can 
arise in condensed matter systems which have 
strong interactions between electrons. Mott insu- 
lators with one electron per unit cell are intrin- 
sically strongly interacting, and a good place to 
look for paradigm shifting phenomena. In particu- 
lar, if no symmetry is broken, the ground-state is 
guaranteed to be an exotic "spin liquid" . We have 
discussed two classes of spin liquids: topological 
and critical. The former exhibit topological order 
and gapped excitations with quantum number 
fractionalization. The latter display power-law 
correlations in the ground state and gapless exci- 
tations. In both cases, gauge theory offers a simple 
way to characterize these spin liquids and their 
excitations, and to encode the statistical interac- 
tions between them. Duality transformations also 
let us gain understanding of these novel quantum 
phases, which in some cases can be probed experi- 
mentally, as in the spin 1/2 compound CS2CUCI4. 

We have also explored in the last section paths 
beyond Landau's theory of phase transitions 
through the example of the deconfined quantum 
critical point between a Neel and a Valence Bond 
Solid state in two dimensions. We reviewed the 
recent theoretical work of Ref. [32] that proposes a 
continuous transition between these phases, with 
emerging deconfined spinon excitations right at 
the critical point. 

These recent developments hopefully convey the 
current excitement in the field of strongly interact- 
ing quantum condensed physics. The specific Z% 
and U(l) spin liquids that were discussed above 
surely represent only the tip of a large iceberg con- 
sisting of many other exotic states of matters that 
can arise, especially in doped systems. But much 
hard work is needed to further understand and 
(hopefully!) detect such phases in the laboratory. 
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